Molecular spin-orbit excitations in the Jofr = 1/2 frustrated spinel GeCo204 
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We describe powder and single-crystal inelastic neutron scattering experiments on a spinel-type 
antiferromagnet GeCo204, represented by an effective total angular momentum Jch = 1/2. Several 
types of non-dispersive short-range magnetic excitations were discovered. The scattering intensity 
maps in Q space are well reproduced by dynamical structure factor analyses using molecular model 
Hamiltonians. The results of analyses strongly suggest that the molecular excitations below Tn 
arise from a hidden molecular-singlet ground state, in which ferromagnetic subunits are antiferro- 
magnetically coupled. The quasielastic excitations above Tjv are interpreted as its precursor. A 
combination of frustration and J^a = 1/2 might induce these quantum phenomena. 

PACS numbers: 75.30.-m, 75.40.Gb, 75.50.Xx, 75.50.-y, 78.70.Nx 



I. INTRODUCTION 

Since their initial proposals,A"" the concept of geomet- 
rical spin frustration has been intensively studied. Geo- 
metrical frustration has been shown to give rise to novel 
forms of spin-liquid-like fluctuations in a paramagnetic 
phase, such as spin molecules, spin ices, and spin vor- 
tices Recently, dynamical spin molecules were discov- 
ered as non-dispersive excitation modes within a mag- 
netically ordered phase, where frustration was assumed 
to be relieved by a lattice deformation.^ 

Meanwhile, it was demonstrated that an effective to- 
tal angular momentum Jeff = 1/2, generated by a spin- 
orbit coupling (SOC), provides a new playground for 
correlated electrons. For example, a Mott instability 
with spin-orbit integrated narrow band was confirmed in 
Sr2lr04, and a quantum spin-hall effect at room temper- 
ature was theoretically predicted for Na2lr03.^°"^^ These 
iridates possess Ir'*"'" with low-spin (t2g)^ configuration, 
of which the ground states are described by Jeff = 1/2 
with unquenched orbital angular momentum {L — 1). 
The L = 1 states are related to t2g triplets {xy,yz, zx): 
= ±1) = {\yz)±i\zx))/V2 and = 0) = \xy). The 
value of 1/2 and the complex orbitals of J^e are expected 
to enhance the quantum nature accompanied with orbital 
degree of freedom. 

Then, an interest in the combination of frustration and 
Jeff = 1/2 will naturally arise. In fact, there are re- 
ports of Ir'*"'" systems, e.g. a face-centered cubic system 
K2lrCl6 forming a magnetic complex IrClg with remark- 
ably mixed orbital, a hyperkagome system Na4lr308 with 
quantum spin liquid, and pyrochlore systems Ln2lr2 07 
(Ln=Nd, Sm, Eu) with metal-insulator transitions^"— 
However, an extremely strong neutron absorption of Ir 



nuclei (~ 425 barns for thermal neutrons)^^ and lack of 
large single crystals hamper the successful inelastic neu- 
tron scattering experiments, a prime tool for the study 
of magnetic frustration. 

The spinel- type antiferromagnet GeCo204 is a promis- 
ing candidate with frustration and Jeff = 1/2. In this 
material, well-known SOC-active Co^+ ions octahedrally 
surrounded by anions form a lattice of corner-sharing 
tetrahedra, which is geometrically frustrated, and Ge''+ 
ions are nonmagnetic. Figure[TIa) shows the energy-level 
schemes of a single-ion state of Co^+ {d"^)- The crys- 
tal field and SOC yield J^e = 1/2, 3/2, and 5/2 states 
with L = 1 and S = 3/2ji2r— Antiferromagnetic order 
with propagation vector q„ = (1/2, 1/2, 1/2) and a tiny 
tetragonal lattice deformation (c/a ~ 1.001) simultane- 
ously occur at Tn ~ 21 K, which is suppressed compared 
to the Curie- Weiss temperature 6w — 81 K.— "— Spin- 
liquid-like fluctuations above Tjv (quasielastic mode) and 
a non-dispersive magnetic excitation mode below T/y (4- 
meV mode) were also found by powder inelastic neutron 
scattering!^ 

In this paper, we comprehensively study magnetic exci- 
tations above and below Tn in GeCo204 in wide momen- 
tum (Q) and energy (E) ranges by powder and single- 
crystal inelastic neutron scattering. The experimental re- 
sults and numerical analyses strongly suggest manifesta- 
tion of highly-frustrated quantum states in this cobaltite. 



II. EXPERIMENTS 

Initial single crystal studies were performed at the 
NIST Center for Neutron Research using the BT-2 and 
BT-9 triple axis spectrometers. Single-crystal inelastic 
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FIG. 1: (Color online) (a) Energy level scheme of Co^^ ion 
under octahedral crystal field and SOC. (b) Correspondence 
with the data measured in the present experiments. 



neutron scattering experiments were performed on the 
triple axis spectrometer TOPAN, installed at the JRR- 
3 reactor, JAEA, Tokai, Japan. The final energy of the 
neutrons was fixed at Ef — 13.5 meV with horizontal col- 
limation sequence of blank-lOO'-IOO'-blank. A sapphire 
filter and a pyrolytic graphite filter efficiently eliminated 
fast neutrons and the higher order contamination, respec- 
tively. Single-crystal rods of GeCo2 04 were grown by a 
floating zone method. Details of the crystal growth are 
summarized in Ref^-. The rod size was about 4 mm di- 
ameter and 30 mm height. The three co-aligncd single 
crystals were enclosed with He exchange gas in an alu- 
minum container, which was placed under the cold head 
of a closed-cycle He refrigerator. 

Powder inelastic neutron scattering experiments were 
performed on the direct geometry chopper spectrome- 
ter HET, installed at the spallation neutron source, ISIS 
Facility, UK. The energy of the incident neutrons was 
flxed at Ei = 59 and 29 meV. A 35 g powder specimen 
of GeCo2 04 was synthesized by a solid state reaction 
method, filled in an envelope made from thin aluminum 
foil, and inserted in a refrigerator with He exchange gas. 



III. RESULTS 

Figures [D^a) and (Hb) show the powder data with 
Ei = 29 meV. Above Tn the quasielastic mode is ob- 
served around Q = \Q\ ~ |q„J ~ 0.66 A~^, as shown 
in Fig. [5Ja). Below Tjv spin- wave-like dispersion rises up 
from around Q — \ \ in addition to the previous discov- 
ered 4-meV mode?^ as shown in Fig. [H^b). Figures HJc) 
and[2jd) show the data with Ei = 59 meV. Two discrete 
levels are discovered around E = IQ and 29 meV both 
above and below T/v, indicating that the two modes are 
not spin waves. Below Tjv these modes slightly sharpen 



and harden. 

We measured Q correlations of the quasielastic mode 
above Tm and the 4-meV, 16-meV, and 29-meV modes 
below Tjv in a constant-i? scan mode by single-crystal 
inelastic neutron scattering, as shown in Figs. EJa) to 
[31[h). The scattering intensity distributions with char- 
acteristic patterns decrease at higher Q, as expected for 
the Co magnetic form factor, indicating that the exci- 
tations must be attributed to magnetic origin and not 
phononic. Figures [3Ja) and ^^h) show the data for the 
quasielastic mode, measured ai E = A meV. The inten- 
sity is strong only in the 400, 440, and 222 Brillouin 
zones, and is distributed near the edges of the zones. 
Figures [U^c) and EJd) show the data measured at 4 meV 
below Tat . Though it is difficult to remove the spin- wave 
component spread around h/2 k/2 1/2 reciprocal lattice 
points (magnetic Bragg reflection points for elastic scat- 
tering) , the scattering pattern is quite similar to that for 
the quasielastic scattering. Figures [Sje) to|3l^h) show the 
data for the 16-meV and 29-mcV modes. The scattering 
intensity of the former mode is relatively strong except 
for the above Brillouin zones, whereas that of the latter 
mode is distributed on every zone boundary. 



IV. MODEL ANALYSES 



We analyzed the quasielastic mode using a molecular 
model, as for spin frustrated systems-^'^ For elastic and 
quasielastic magnetic neutron scattering, the cross sec- 
tion is described by 



N 



S{Q) - Co|F(Q)|2|^ J,x exp(zg • r,)\ 



(1) 



where Co is a proportional constant of intensity, F{Q) 
is the magnetic form factor of the Co^+ ion, for which 
the Watson-Freeman one was used below j labels the 
site of the Co^"*", N is total number of the sites in a 
molecule, rj and rji are those positions, and Jj± is an 
expected value of Q-perpendicular component in Jj-'^ 
When colinear J's fluctuate in arbitrary directions like 
in a hexagonal-type quasielastic mode observed in the 
typical spin-frustrated system ZnCr204, Jj± (= Sj±) 
takes on only ±1.^ Following this treatment, we searched 
for and found a di-tetrahedral model for the quasielastic 
mode in CeCo204, as shown in Fig. UJ^a). Figures [3l^i) 
and[3Jj) show the calculated patterns, which are in good 
agreement with the experimental patterns of Figs. [3Ua) 
andEIb). 

For inelastic magnetic neutron scattering, the cross 
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FIG. 2: (Color online) (a)-(d) Color images of powder inelastic neutron scattering data with different incident energies and 
temperatures. The color tones indicate the scattering intensity in mbarn/(sr-formula) units, (e) Energy spectra, averaged from 
3 to 29 deg in scattering angle (Q — 0.3 to 2.7 for elastic condition) in (c) and (d). The arrows indicate the 16-meV and 
29-meV modes. 



section is described by 

co\Fm's{f^ - ^) ( E ('^"^ - X 

N \ 

^ (A I I A')(A' I J| I A) exp{zQ • (r, - r,,)} , 

(2) 

where a and /? are {x,y,z), N is number of sites in a 
molecule, | A) and | A') are molecular ground and excited 
states, respectively, J is a total angular momentum op- 
erator, and the parentheses indicates an orientational av- 
erage over equivalent molecules.— One cannot generally 
apply Eq. (IlJ, which is obtained from Eq. ([2]) only when 
the matrix elements (A' | J" \ A) can be simply reduced 
for elastic scattering. In the following we try to reproduce 
our inelastic scattering data using a relatively simple 
molecular model. We assume effective molecular Hamil- 
tonians, and numerically evaluate (A' | J" | A) and the 
cross section. The assumption of a molecular formation 
implies a remarkably-mixed molecular orbital, which will 
enhance intra-molecular exchange interactions and sup- 
press atomic orbital characters like anisotropy^SS There- 
fore, we ignore the exchange field outside the molecule 
(Lorentz-like local magnetic field) and the directional 
term {6a js — QaQp/lQl"^) in Eq. ([2]). An orientational 
average over dynamically fluctuating molecules will also 
substantially suppress the directional dependence. For 
simplicity the atomic Watson-Freeman form factor is 
used for F{Q) in Eq. © againM 

Firstly, we exactly diagonalized a tetramer Hamilto- 
nian: 

4 

-fftctra — Ji ' Ji ■ J j, (3) 



where Ji = 1/2, i and j are positions of the tetra- 
hedral sites (Fig. IHJd)), J2i j means summation over 

all J pairs (not doubly counted), and is a first- 

neighbor exchange interaction that is ferromagnetic as 
expected from the quasielastic mode (Fig. IDJa)) and the 
Goodenough-Kanamori rule.^° The 16 (= 2^) basis states 
of \Jl, J2, J§, Ji) were used, where Jf = ±1/2. Fig- 
ureSl^d) shows the obtained level scheme with J^^^"^ = —8 
meV. The ground states are described as ferromagnetic 
quintets with Jtctra — 2, and the first excited states 
are nonets with Jtetra = 1 and E = 16 meV, where 
•Jtctra — J2i=i '^i- The nonet can generate all states with 
"^tetra = ±1,0 with a Jeff = 1/2 dimer-singlet bond by 
their linear combinations (e.g. Fig.|4l^d)). The calculated 
patterns for excitation processes from the ground states 
to the excited states are shown in Figs. [3l^m) andjU^n), 
which are in excellent agreement with the experimental 
patterns ofFigs.EJe) and[3Jf) (16-meV). 

Secondly, we diagonalized a di-tetramer Hamiltonian: 

8 

i?di-tetra = j\f ^ ^ ' ' (4) 

where J^^^^ — j'f^'^ and J^^\ i and j are positions of 
sites in the di-tetramer (Fig. SJc)), and the 256 (= 2*) 
basis states of | Jf , J|, J|, J|) were used. The sign 
of J^'^^'^ is antiferromagnetic, being consistent with the 
quasielastic mode and previous neutron diffraction re- 
ports. Figure m^c) shows the level scheme with J^'^^^ — 
—8 meV and jj^'^^'' = 10 meV. The ground state is de- 
scribed as a non-magnetic singlet with Jdi-totra = 0, and 
the first excited states are triplet with Jdi-tctra — 1, where 
Jdi-tctra = J2t=i Ji Figures EJk) andlSJl) show the cal- 
culated patterns of the singlet-triplet excitations, which 
are similar to those for the quasielastic mode (Figs. [Sji) 
and[3Ij)), and are identified as the 4-meV mode. 

For the 29-meV mode, we could find no model within 
Ji = 1/2 after many trials. On the other hand, in- 
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FIG. 3: (Color online) (a)-(li) Color images of single-crystal inelastic neutron scattering data, measured in the hkO and hhl 
zones in a constant energy scan mode, (a) and (b) were measured at iJ = 4 meV. (i)-(p) One-to-one correspondence between 
calculated patterns as identified by the molecular models shown in Figs. [U^a) to and described in the text. The bold 
lines show the Brillouin zone boundary of the spinel structure. For the calculated patterns, the horizontal bars indicate the 
scattering intensity in arbitrary units. 



terestingly, other cobalt compounds KC0F3, CoO, and 
La(Sr)2Co04, consisting of Co^+ ions octahedrally sur- 
rounded by anions as well, exhibit excitations around 
30 meV^i^i^ These excitations are interpreted as the 
lowest-energy SOC excitations {i.e. excitons)i^"— In 
analogy with these cobalt compounds, the 29-meV mode 
in GeCo204 is to be excitons. 



Thus, thirdly, we studied molecular excitons from 
^off = 1/2 to Joff = 3/2, assuming the following Hamil- 
tonian of third- neighbor dimer (Fig. UJe)): 

H,,^jt^'j2s.-S„ (5) 
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FIG. 4: (Color online) (a)(c)(d)(e) Schematic representations of the molecular models. The green arrows represent magnetic 
moments of the Co'^"'' ion, and green ellipsoids represent a non-magnetic singlet formation. All the moments dynamically 
fluctuate in arbitrary directions with the relative correlations. The structural units shown in (a) and (c) are identical to each 
other. In (c) to (e), the representative states are depicted, (b) First and third neighbor exchange interactions. Representative 
bonds are shown. 



where j!^'^^'' — 2 meV is an effective value of j!f^'' = 10 

meV for spin-3/2 estimated by a relation J3'''''''5(5+l) = 

J^'^^^ Jcs{JcS + 1), and i and j label sites in this dinier. 
This Hamiltonian is expressed by S, not J, for the 
cobalt compounds, The 36 {— 6^) basis states of 
I Ji, Jf ) ® I J2, J|) were used, where 1 1/2, ±1/2) single-ion 
states have zero energy, and |3/2, ±3/2) and |3/2,±l/2) 
ones have a 29 meV SOC excitation energy for \ Ji, Jf). 
Figures[31[o) andj^Jp) show the calculated patterns, which 
take into account the processes from the Jeff = 1/2 
ground dimer-singlet to the first excited triplet with 
JcS = 3/2. The calculated patterns are in good agree- 
ment with the experimental patterns (Figs. El^g) and 
Hh)). 

More precisely, the Co^+ feels an additional trigo- 
nal component of crystal electric field, which keeps the 
JcS = 1/2 ground doublet but splits the Jeff 3/2 
quartet into two Kramers doublets (J^jj = ±3/2 dou- 
blet and Jgff = ±1/2 one) in the level scheme shown 
in Fig. [IJa)i ^^i^° In fact, as shown in Fig. [DJe), the ex- 
perimental spectrum around 29 meV is asymmetrically 
spread up to 40 meV, suggesting this splitting. However, 
the profile is too broad to clearly resolve into the two 
levels. Therefore we carried out the above calculation in- 



tegrating the Jeff — 3/2 quartet. We also confirmed that 
the transitions to each Kramers doublets give the same 
patterns. 

In this way, we identified the quasielastic mode as 
an antiferromagnetic di-tetrahedral cluster (Fig. HJ^a)), 
consisting of Co^+ ions with Jeff = 1/2. Further- 
more, assuming model Hamiltonians, we assigned the 4- 
meV to the singlet-triplet excitations in a di-tetramer 
with the same structural unit (Fig. IDJc)), the 16-meV 
to quintet-nonet excitations in the one ferromagnetic 
tetramer (Fig. SJ^d)), and the 29-meV to SOC excitons 
from Jeff = 1/2 to 3/2 in a third- neighbor-distant antifer- 
romagnetic dimer (Fig. |4je)). All the excitations can be 
regarded as intra-activations of the di-tetramer. The cor- 
respondence relation between excitations and Jefi states 
are shown in Fig. W[h). 



V. DISCUSSION 

We discuss the ferromagnetic tetramer (Fig. IH^d)). 
The remarkable spatial confinement of magnetic corre- 
lation demonstrates the existence of frustration. But 
frustration is normally based on antiferromagnetism. So 
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what is frustrated in GeCo204? One factor will be the 
frustration among J^'^^^ , jj^'^^^ , and the other exchange in- 
teractions. Aside from this, however, an orbital system is 
inherently frustrated even on a simple cubic lattice; when 
orbitals (directions of the electron cloud) are arranged to 
gain bond energy for one direction, this configuration is 
not fully favorable for other bonds,— GeCo204 also 
has an orbital angular momentum, which is a kind of or- 
bital, and is in the geometrically frustrated pyrochlore 
lattice. Therefore, both exchange and orbital frustration 
likely coexist in GeCo204. 

Next we discuss the di-tetramer. According to the 
above analyses, the di-tetramer singlet ground state is 
surprisingly hidden as origin of the molecular excitations 
below T/v- Indeed, a singlet formation is an effective 
way to suppress frustration and degree of freedom. How- 
ever, the formation does not necessarily mean that all 
the magnetic moment disappears, because the g factor is 
arbitrary in our analyses, being consistent with the coex- 
istence of singlet and magnetic order. This partial-singlet 
model can explain why GeCo204 exhibits magnetic order 
with only about 3 /ie per Co^+,'^^ which is 1 fiB lower 
than a normal value of 4 fiB generated by SOC like in 
CoO^ 

It should be noted that a typical spin-frustrated spinel 
antiferromagnet MgCr204 (Cr^+, (P, S — 3/2) similarly 
exhibits a set of quasielastic modes above Tjv (hexamer) 
and a gapped non-dispersive excitation mode below T/v, 
of which the scattering intensity distributions in Q space 
are the same.- In addition, MgCr204 exhibits magnetic 
order with only 2.2 /iB,"^^ which is about 1 //b lower than 
the full moment 3 fiB- Therefore, a hexamer- type singlet 
ground state would give rise to both the gapped excita- 
tion mode and the partial disappearance of the magnetic 
moment below T/y. We also remark that the 1 /iB de- 
crease is observed in the isomorphic systems ZnCr2 04 
and }lgCT204^ 

The di-tetramer can be energetically regarded as a 
dimer of the rigid tetramers with Jtetra = 2; bind- 
ing energy in a ferromagnetic tetramer (~ 36 meV = 
gj(ex) ^ Jgjj^j^jj -)- 1)) is higher than antiferromagnetic 
coupling energy between the two tetramers (~ 24 meV 

= (44°''^ + J^"""^) ■ Jeff(Jeff + 1)). We also numerically 
confirmed that the Jtetra- 2 dimer has a ground singlet 
with the combination of Jtcu-ai — ±2, ±1,0 and the 
first excited triplets within the 25 {— 5^) basis states 
of I Jtetra,!' •^tetra,2)- This cxtendcd-dimer picture natu- 
rally gives us the interpretations of the 4-meV mode as 
a localized singlet-triplet excitation and of the quasielas- 
tic mode as its precursor fiuctuations, as observed in the 
frustrated spin-1/2 system SrCu2(B03)2 with the two- 
dimensional Shastry-Sutherland latticci^ 

Since its introduction as a mechanism for high- 
temperature superconductivity, dimer-based quantum 
cooperative phenomena like resonating valence bond 
(RVB) and valence bond solid (VBS) have been sought 
after in fields of magnetism and strongly correlated elec- 



tron systems.^ SrCu2 (303)2 is one of the great suc- 
cesses. In contrast to the borate, the molecular for- 
mations in GeCo204 are characterized by the existence 
of a ferromagnetic molecule, Jes = 1/2, and the three- 
dimensional pyrochlore lattice with almost regular tri- 
angles. In this sense, GeCo2 04 could be positioned as 
a new class of quantum cooperative systems caused by 
frustration and JoS = 1/2. 

We list two other intriguing characters of frustration 
and Jeff = 1/2. One character will be the fact that all the 
molecular excitations involve aspects of not only spin but 
also orbital excitations by SOC (molecular orbitons). An- 
other character will be the emergence of molecular exci- 
tons (29-meV). Excitons normally appear within a single 
atom or ion with SOC, and are occasionally propagated 
with very narrow dispersion width (^0.5 mcV) like in a 
4/ electron system4i Furthermore, the 3d electron cobalt 
systems exhibit more dispersive excitons around 30 meV 
(over 5 meV width), propagated by stronger exchange 
interactions than in 4/ systems<^"— However, these ex- 
citons are molecular, which are locally collective but are 
not propagated. 



VI. CONCLUSIONS 

We discovered several types of non-dispersive short- 
range excitations in a three-dimensional frustrated 
GeCo204 with J^s = 1/2 by powder and single-crystal 
inelastic neutron scattering. The scattering intensity 
maps in Q space are well reproduced by quantum- 
mechanical molecular models. The model analyses 
strongly suggest that a molecular-singlet ground state 
consisting of ferromagnetic sub-molecules is hidden be- 
low T/v, which gives origin to the molecular excitations. 
The quasielastic excitations above T/v are interpreted as 
a precursor of this quantum ground state. The spin and 
orbital frustrations of Jcs lead to the molecular-singlet 
formation and the ferromagnetic molecule one, respec- 
tively. Further experimental and theoretical works will 
be needed to fully elucidate this hidden molecular partial- 
singlet conjecture and clarify the molecular orbital for- 
mations. 
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